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We study experimentally the influence of the intensity of the infrared (IR) probe field on attosec- 
ond pulse train (APT) phase measurements performed with the RABITT method (Reconstruction 
of Attosecond Beating by Interference in Two-Photon Transitions). We find that if a strong IR field 
is applied, the attosecond pulses will appear to have lower-than-actual chirp rates. We also observe 
the onset of the streaking regime in the breakdown of the weak-field RABITT conditions. We per- 
form a Fourier-analysis of harmonic and sideband continuum states and show that the mutual phase 
relation of the harmonics can be extracted from higher Fourier components. 
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FIG. 1: Setup for the experiment: a 2mJ, 30 fs laser pulse is split into a probe and a pump arm. The larger part is used in 
high-order harmonic generation while a fraction travels through a delay stage to serve as probe pulse. The generated harmonics 
are filtered spectrally and spatially by silica plates, a thin-film metallic filter and an aperture, before recombining with the 
probe using a drilled mirror and being focused into the detection gas of a magnetic bottle electron spectrometer. 

I. INTRODUCTION 

Laser-assisted ionization processes provide an elegant tool to study dynamics and details of atomic and molecular 
systems [2[5J[3]. For time-dependent measurements, the assisting laser needs to be synchronous to the ionization 
event on a time scale shorter than that of the process to be resolved. A number of experiments in the past years have 
shown that the attosecond (as) time scale is routinely accessible and laser- assisted ionization processes have become 
the chief tool for experiments in the field of attosecond science [U El [3l HJ |5j |6] . 

The emission of a comb of high-order harmonics when an atomic medium is exposed to a driving intense laser field 
is well understood. The resulting attosecond pulse trains (APT) |S] provide a premier tool to controllably ionize 
atomic media [9] and it is important to determine their characteristics, both amplitude and phase. This can be 
done by cross-correlation with an IR probing field under stable interferometric conditions. Using weak infrared (IR) 
fields, this is often done with the Reconstruction of Attosecond Beating by Interference of Two-Photon Transitions 
(RABITT) method [7 , while the AC-Streak camera method [10, 11., using stronger IR fields, has mainly been applied 
to the characterization of isolated attosecond pulses. A special case is the full characterization method FROG-CRAB 
(Frequency- Resolved Optical Gating- Complete Reconstruction of Attosecond Bursts) , which is based on an iterative 
deconvolution of a time-frequency spectrogram |12l 113] . In this paper, we study the transition from the weak to the 
strong field regime, and the effect of the probe field strength on the validity of the characterization. 

In Section II, the experimental setup and performed experiments will be presented. We perform in Section III 
a conventional analysis of the data in the RABITT regime, studying specifically the probe intensity dependence of 
the phase measurement. A Fourier-series approach (Section IV) allows us to see the fingerprints of processes with 
more than one contributing IR photons and provides a more general description of our delay-dependent laser-assisted 
photoionization spectra. 

II. EXPERIMENT 
A. Setup 

The experiments were performed at the High Power Laser Facility of the Lund Laser Centre. We use a 1kHz 
chirped-pulse-amplification (CPA) titanium-sapphire based laser system, providing 2mJ 30 fs pulses centered around 
800 nm. The pulses are focused into an Ar gas cell, pulsed at 1 kHz, to generate high-order harmonics, which then 
propagate through silica plates and a thin Al-foil to filter the harmonics to form a well-defined APT (cf. Fig. [l}. 
The APT is then recombined with an IR probe beam in a Mach-Zchnder-type interferometer. In the last step of the 
interferometer, the APT is passed through a drilled mirror while the IR is reflected on the back side of the mirror. 
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FIG. 2: XUV photoionization spectra, recorded for varying sub-cycle delay between IR and XUV field. The corresponding 
probe pulse energies are written in the upper part of the figures. A change of this energy translates proportionally to a change 
of the probe intensity. The zero delay indicates the delay at which the maxima of the IR pulse and the APT coincide. 



Both beams are focused into the detection gas of a magnetic-bottle electron spectrometer (MBES) with the help of 
a toroidal mirror. Using a combination of half- wave plate and polarizer in the probe beam, the pulse energy of the 
IR can be varied continuously between 5 and 100 /iJ. As all focusing parameters are kept constant, any change in 
the probe pulse energy will directly result in a proportional change of the IR intensity in the detection region of the 
MBES. The MBES has a 2tt acceptance angle and a maximum energy resolution of about 100 meV. 

The relative delay of IR probe and APT can be adjusted on two time scales. Using a motorized translation stage 
in the probe arm we can vary the delay in the range of one femtosecond to several picoseconds. To accurately resolve 
attosecond processes, a delay stage with a piezoelectric crystal is used to change the relative phase of the two beams 
with a precision of a few tens of as. This stability allows us to perform interferometry of two-photon pathways as in 
the RABITT characterization scheme, as shown below. It is also crucial for the use of the AC streak camera method, 
where the relative phase of the two beams has to be stable with similar precision. 



B. Delay-Dependent Two-Color Photoionization Spectra 



In the experiments presented here, the resulting photoelectron spectra from an IR probe and XUV field were 
recorded. A shift in the relative delay of the two fields induces a change in the observed spectrum, and recording the 
spectra at different subcycle delays allows us, in some conditions, to determine the temporal structure of the XUV 
emission. 

Figure [2] shows the seven scans that constitute our experimental results. From left to right, the probe pulse energy 
changes by a factor of 20, ranging from 5 to 100 fi J. At the lowest intensity, the presence of a weak IR field leads to 
weak sidebands between the odd harmonic orders. These sidebands originate from u>„ + uir and ui n — lor two-photon 
transitions [2] where lor is the IR photon frequency and u>„ = nojji, with n being an odd integer, see Fig. [3] a). The 
two possible pathways to each final state lead to the observed interference pattern. This pattern is repeated every 
IR half-cycle, Tr/2, reproducing the frequency of the attosecond pulse periodicity. This probe intensity regime is 
commonly referred to as the RABITT regime. 

With increasing IR intensity, the amplitude of the sidebands becomes comparable to direct photoionization by the 
harmonics [see, e.g. sideband 22 in Figure [2ji)-g)]. These two-photon processes induce a depletion of the peaks at odd 
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FIG. 3: Sketch over possible transitions involving single harmonic- and multiple IR photons, a) shows the RABITT transition, 
which couples two harmonic photons with an energy difference of 2ljr to a final state with energy (n — 1)ljr. b) shows that the 
absorption/emission of two IR photons leads to a similar coupling with odd final energy tilor. c) illustrates that yet another 
IR photon can be used to couple harmonic photons that are 6u>r apart in an even final energy (n — 1)u)r. d) and e) shows 
that three IR photons will also couple new harmonics that are 2lor apart in even final energy (n — 1)wr, which will distort the 
information from the RABITT signal in a). 



harmonic energies, at the delays where the sidebands are maximum. The cutoff region is first affected by depletion 
since the strength of continuum-continuum transitions rapidly increases with the energy of the initial state [Fig. 
and c) at harmonic 25 and above]. For a probe pulse energy of 40 /zJ [Fig. |2ji)], depletion effects become visible even 
in the low energy region [see harmonic 19 to 23 in Fig. [2ji)-g)]. 

At the higher probe intensities [Figure [2ji)-g)], processes involving more than one IR photon become significant. 
This is the so-called streaking regime, where the AC-streak camera becomes the preferred characterization method 
PHI mi H^] . The streaking regime is clearly entered in the last two of the scans in Fig. [2j the cutoff being increased 
by as much as seven IR photons (from harmonic 27 and above) . Note that the electron signal is still showing discrete 
IR photon energy spacing because of the periodicity of photoionization from the sequence of attosecond pulses. 



III. ANALYSIS OF THE EXPERIMENT 



A. Reconstruction of Attosecond Beating by Interference in Two-Photon Transitions (RABITT) 



We now analyze the results in Figure [4] [same as|2p)] within a perturbation theory framework. Here, we will consider 
the photoionization of the detection gas by the APT as a first-order perturbation and the onset of the sidebands and 
the RABITT modulation as a perturbation of the second order. This is a good description of two-color photoionization 
at low IR intensities. 

Contributions to the final state probability amplitudes of order N can be obtained from those of order (N-l) through 
the following equation [TBI H6j : 



5 



32 



30 



0) 

"D 
i_ 

o 

o 
'c 
o 

E 
i_ 
ro 
X 



28 
26 
24 
22 
20 
18 
16 



14 



Delay (fs) 



FIG. 4: Delay-dependent photoionization spectrum in the RABITT regime, 
line indicates the sideband maxima. 
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a ^'(t) = ^EE / dt'^^M^it'v^ 1 -^', (i) 

where I is the index of the respective contributing states, u>f are the frequencies acting on the atom and m denotes the 
final state. The effect of the IR on the ground state is negligible. The population of continuum states after perturbing 
the system with a harmonic frequency can be aproximated by 

a«(t) = Vn, 9 E(u Jn ) t (2) 

where fi n g is the transition dipole matrix element from the ground to the continuum state, and t the interaction time, 
approximately a measure of the pulse duration. The population of the sideband state (n — 1) can be expressed as 

1 t 2 

fl L lW = J^2^ [Vn-l,n-2^n-2,gE{uJn)E(uJ n _2) + Hn-l,nHn,gE{-UJ R )E(uJ n )\ (3) 
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(ih) 2 2 



where the factor M represents a combination of the various transition dipole moments, which we approximate as 
equal for both quantum paths, and E R0 , E X o ar e the amplitudes of the IR and XUV fields. The sideband intensity 
will finally oscillate as 

&(<p) « cos (24> R + A^ 2i „), (4) 

with (f)R — lurt, t being the delay between the XUV and the IR fields. A0„_ 2 ,n is the phase difference between 
the harmonic spectral phases <fi n -2 and <f> n . This allows us to obtain the phase difference A^„_2 ! n from a Fourier 
transform of the spectrum over delay. 

Using the obtained phase, we can reconstruct a pulse shape as in Figure [5j The Fourier limit for our spectrum 
is 160 as. Due to the chirp rate of 18400 as 2 we obtain an average pulse duration of 440 as FWHM. This method 
of measuring APTs has been successfully employed in a number of experiments [7J [51 [T7] . The high chirp rate and 
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FIG. 5: Reconstruction of average pulse in the pulse train, measured by the RABITT method. The upper panel shows 
the harmonic intensities (solid line) and integrated spectral phase (dashed line). Using these intensities and phase, one can 
reconstruct the pulse shape as in the lower panel. The Fourier limited pulse shape is shown as dashed. 

asymmetric pulses come as no surprise as only one thin aluminum filter was employed to counter the intrinsic chirp 
of the HHG process. The addition of more filters would help to approach the Fourier limit and further compress the 
pulses [T7J [TB] . The RABITT method provides access to the relative phase difference of presumably monochromatic 
harmonics [TJ [14] . The measured phase difference is equal to the group delay (GD) over the spectrum of the APT, and 
an integration allows us to reconstruct an average attosecond pulse in the train. Such a retrieved pulse represents a 
good approximation for pulses within the FWHM of the APT, accounting for about 90% of the signal from experiments 
with such trains. There are extensions and implementations of the method that yield more information on the full 
structure of the APT [191.120]. 

B. Intensity Dependence of RABITT Signal 

The influence of an intense probe field on the measured spectral phase was studied by progressively increasing 
the IR intensity and recording the delay-dependent photoionization spectra. Figure [6] shows the phase differences 
A</>„_2,n evaluated according to the method outlined above, for varying probe intensities. It is clear that the slope is 
decreasing with increasing intensity, as if the chirp rate was lower, even though the actual APT remained the same. 

As a consequence of the phase variation with intensity shown in Figure [6] the reconstructed pulses will appear 
more compressed at higher probe intensities, which is a pure artifact of the increased IR intensity. The RABITT 
analysis relies heavily on the presence of a single IR photon contributing to the final state as to allow only the phase 
difference of the neighboring harmonics to be measured. Higher-order processes perturb the RABITT signal, as will 
be discussed in more detail in Section ITVl 

C. Streaking 

We now consider briefly the regime of high probing intensities [Figure[7] same as[2]g)] where the interaction between 
photoelectrons and probing laser field can be understood by simple classical arguments. In this regime, an electron is 
released into the continuum at a given ionization time ti due to photoionization by an attosecond pulse. It will gain an 
additional momentum proportional to the IR vector potential A(ti) at the time of ionization. This momentum shift 
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FIG. 6: Phase differences A(f> measured over the harmonic spectrum, for different probe pulse energies. Higher probe intensities 
result in a flatter measured spectral phase in the plateau region (order 14-22). 



imparted by the IR probe is therefore dependent on the relative delay between the two fields. With APTs generated 
by two-color fields, or with single attosecond pulses, the pulse properties (duration and chirp) can be extracted from 
analyzing the streaking trace and the pulse(s) can be reconstructed [lOj [TTJ [13]. In the case of an APT with two 
pulses per cycle, the situation is more complex and a reconstruction of the attosecond pulses from the experimental 
data is more difficult, since the sign of the probe field reverses for consecutive pulses. Without going further with 
the analysis of our experimental data, we note that the chirp of our attosecond pulses, however, is clearly visible in 
the high energy region of the spectrum [see dashed line in Figure [7] . The onset of the streaking regime also explains 
the observed phase shift of the sideband orders close to the cutoff in Figure [6] In the perturbative regime and, for 
simplicity, in absence of attosecond chirp, the maxima of the sidebands occur at the maxima of the probe laser electric 
field. At moderate probe intensity, depletion of the harmonic states will lead to maxima of the harmonics shifted by 
7r/2 with respect to the sidebands, thus corresponding to the zeros of the electric field. At high energy, electron peaks, 
coming from direct electrons being streaked by the probe field, will therefore be maxima at the zeros of the electric 
field, thus shifted by a factor 7r/2 from the sidebands. 



IV. FOURIER DECOMPOSITION OF THE PHOTOELECTRON SPECTRA 

A. Perturbation Theory Analyzis 

To better understand the influence of a high IR intensity on a RABITT measurement, we extend the analyzis 
performed in the previous section, based on perturbation theory to the next orders. The third order of perturbation 
includes components coupling the states {n + 2) and (n — 2) to the harmonic (n) state by two IR photon absorption 
or emission [Figure pH: 
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FIG. 7: Delay-dependent photoionization spectrum in the streaking regime, measured at 100 /^J probe pulse energy. 



This term gives rise to oscillations at frequencies 4ujr, 2u>r, and to a constant term. The intensity of the Aljr 
modulation at the harmonic energy ulor is 

n— 2,n+2 ) (6) 

which is similar to Equation [4j The next order of perturbation leads to a Qujr component in the sideband states and 
the next to an 8ajR component in the harmonic states. Higher IR intensities lead to new couplings of states lying 
further apart and thus higher modulation frequencies. 

Figure [3] illustrates how the number of pathways that lead to a given modulation frequency is unique only under 
certain restrictions. Even numbers of contributing IR photons couple to harmonic states while odd numbers couple 
to sideband states. The frequency is meaningful only in the harmonic final states, as shown in Eq. [5] and the 
phase difference that can be determined is that between harmonic states (n — 2) and (n + 2). A Fourier transform 
that isolates this component will allow us to access to the spectral phase in the same way as with the conventional 
RABITT, looking at the 2u>r modulation in the sideband states. Experimentally, due to the resolution of our scans, 
the 6ojr component is the highest resolvable. It is already close to the Nyqvist limit, with only three to four data 
points per period. 

Fig. [HJshows a comparison of the group delay obtained over the spectrum for three different modulation frequencies. 
The usual RABITT is based on the 5 /iJ-scan, the phases of the 4wfl-modulation were obtained at 20 /iJ, and the 
6w#-modulation was present for the central harmonics in the 80 /iJ-scan. Excellent agreement was found between the 
chirp rates obtained with the different methods. Our results are also consistent with previous experiments [5J 117] . 
From the 2w/}-measurement we obtain a chirp rate of 18400 as 2 or - in terms of a group delay - an emission time 
difference At e — A(/)/(2lur) = 87 as. This is very similar to 16800 as 2 (At e — 79 as) obtained from the 4w7?-componcnt 
in the harmonics. The slope of the phase difference in the 6wfl-component is also in good agreement with 16900 as 2 
chirp rate or a At e of 80 as. 



B. Generalization 



Because of the analogy between the conventional RABITT and the possibility to extract the same information from 
higher components we will now introduce a more general expression for the delay dependence of our photoelectron 
spectrum. This expression will hold beyond the breakdown of the perturbation picture. The amplitude of any 
photoelectron peak in the spectrum is periodic with delay (/jr between the IR and XUV fields. It can therefore be 
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FIG. 8: Comparaison of the obtained phase differences for three different frequency components present in the experimental 
electron signal. The conventional RABITT includes contribution from sidebands 14 to 26. The 4cj_R-component has been 
extracted from harmonics 15 to 23 and the 6w_R-modulation was obtained from sidebands 18 to 24. The curves have been 
shifted for better comparison. 



expressed as a discrete Fourier series: 

oo oo 

s(<j> R ,i R )= s fe (w 0i?2fe = ]rs fc (^,/fl) 



k— — oo 



fe=0 



with St = S-k- Considering only the k — ±1 contributions we have 

= 2 5i(Jk) cos(2^ R + arg(S 1 (I R )) 



(7) 



(8) 



which is the common RABITT case, for low IR intensity, with arg(Si(In)) — A0 n _i. For higher IR intensities, higher 
order terms contribute, the argument will change and higher order terms of the sum Sk(lR) will grow revealing higher 
contributing frequencies. In the case of k = 2, 



S 2 {4>r, Ir) = 2 \§ 2 (Ir) J cos (4^ fl + arg(S 2 (I R )) 



(9) 



which can readily be identified as the Aojr component with an argument of A(f> n —2, n +2 for a suitable probe intensity. 
Figures |3jl)-e) show the effect of the participation of three IR photons, and how two non-surrounding harmonics 
contribute their phase difference to sideband n — 1, destroying the unique source of the 2ujr modulation. Thus 
arg(Si(Ift)) changes and the phase observed in the RABITT analysis varies with IR intensity. 

At the same order of perturbation, a 6lor modulation can be constructed with three IR photons coupling harmonic 
n — 4 and n + 2 to sideband n — 1 [cf. Fig. [3]d)]. For sufficiently high intensity, the number of contributing IR 
photons is 3> 1, resulting in a great number of frequencies added to form the delay-dependent signal in the sideband. 
More and more harmonics contribute to the phase, which becomes flat. The oscillation does not uniquely depend on 
the adjacent orders anymore - the requirement for RABITT. This is the streaking regime, where the classical limit is 
reached. 

As we previously studied the effect of the IR intensity on the phase measurement in the conventional RABITT 
method, we perform a similar analysis for our measurements with the 4wj? component. Figure [9] shows the intensity 
dependence of the harmonic phase measured in this component. We find that at 20 to 60 /j.3 the phase is meaningful 
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Harmonic Order 

FIG. 9: Intensity dependence of the measured phase from the 4uir component in the harmonics. The measured chirp rate over 
the spectrum shows good agreement at 40 fiJ with the 2wi{-measurement from the sidebands. It breaks down at higher probe 
energy. 

compared with the conventional measurement. The measurement at 80 [iJ is showing a similar characteristic flattening 
of the phase as the conventional RABITT. The Aojr breakdown occurs for much higher probe intensities. In this case 
the number of contributing photons becomes so high that arg(S2(In)) is affected by ambiguous coupling of states. 
The analyzis of this 4ur component provides an additional check on the measured phase from conventional RABITT 
method. This allows us to assess whether the probe intensities in any experiment have been too high by comparing 
the 2ujr with the Aojr components and possibly even higher orders. 

V. CONCLUSION 

We have studied the influence of IR probe intensity in two-color XUV photoionization experiments. In the laser- 
assisted ionization process that is at the heart of the RABITT method for characterizing APTs, the probe intensity 
needs to be maintained at levels which only perturb the process very slightly. A high probe intensity will greatly alter 
the measured phase relation of the individual high-order harmonics, making it appear flat over the spectrum, as if 
the pulses were compressed. RABITT is only valid in the limit of weak IR intensity. 

In a next step, the increasing number of contributing IR photons was tracked down by Fourier-analysis of sideband 
and harmonic states. As an increasing number of photons allows to couple states further and further away, higher 
modulation frequencies occur. The 4cur component of the harmonic states allows us to obtain their mutual phase 
relation. Also in this case an increasing IR intensity started to affect the phase measurement making the pulses 
appear artificially compressed. We believe these experiments will allow the scientific community to gain a better 
understanding of the IR intensity dependence of laser-assisted ionization of a gas using high-order harmonics. 

This paper is dedicated to the memory of N. B. Dclonc, who was a father of this exciting field of research "atoms in 
strong laser field", with the first "multiphoton ionization" (1965) and "multielectron multiphoton ionization" (1975) 
experimental results. This research was supported by the Marie Curie Intra-European Fellowship (Attotcch), the 
Marie Curie Early Stage Training Site (MAXLAS), the Knut and Alice Wallenberg Foundation, the Swedish Research 
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Council, and the National Science Foundation (Grant No. PHY-0701372). 
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